Change of the Chern number at band crossings 



in 

On 
On 



J. Bellissard* 
Universite Paul Sabatier, Toulouse, France > 



Abstract 

Let H (e, x) be a finite dimensional hermitian matrix depending on the variable x taking 
its values on a 2D manifold and changing with e. If at e = two bands are touching, we 
" give a formula for the change of Chern number of these bands as e passes through zero. 



> 
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2 1 Introduction 

In many problems of quantum physics, one has to compute the band spectrum of some Hamil- 
tonian. This is what happens in solid state physics for the electronic band spectrum. This 
is also the case in molecular or nuclear physics with rotation bands in the high spin limit 
Jl|, 0. In most of these cases, the Hamiltonian describing the dynamics can be reduced to a 
N- dimensional hermitian matrix- valued continuous function x G M 1— > H(x), where M is some 
manifold of parameters. For the rotation bands, M is the 2D sphere, whereas for electrons 
in a perfect crystal, M is a torus. Diagonalizing H(x) for all x's gives rise to x-dependent 
eigenvalues Ei(x), . . . ,E N (x) called the band functions. Moreover, each of the corresponding 
eigenvectors define a line in depending continuously on x, namely a line bundle. One im- 
^ ■ portant question to answer is whether this bundle is trivial or not, namely if one can choose 
the eigenvectors for each x in such a way as to define a univalued function over M. If H(x) 
is real valued this is always possible. However if H(x) cannot be represented in some basis by 
a real valued matrix depending continuously on x, this is not possible in general. In physics 
such a situation happens whenever the dynamics is not time- reversal invariant, for instance if 
some magnetic field is present or if there are spin couplings. The obstruction is measured by 
the Chern classes of the line bundle. Such Chern classes have been related to the Quantum 
Hall effect |3], ^, || for Bloch electrons in a uniform magnetic field. For rotation bands, the 
occurrence of Chern numbers gives rise to the quantization of a mechanical response to a torque 
in molecules or nuclei . In the present work, we will restrict ourself to the case for which the 
manifold M is two-dimensional, in view of its application to physics. 

We assume now that the Hamiltonian depends upon an additional real parameter e varying 
in a neighbourhood of zero, say [—1, +1]. Since a famous result of Wigner and von Neumann 
J7J, it is known that, if H depends smoothly on (x,e), the set of points (e,x) for which two 
bands of H (e, x) are touching has codimension one generically. 
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From now on, we will assume that M is an oriented smooth surface, namely it has dimension 
2, and that we are in this generic situation. This means that H depends smoothly on (x,e), 
and that there is a pair E-(e, x) < E + (e, x) of neighbouring bands and a finite set Xx, . . . , x^ of 
points in M such that E_(e, x) = E + (e, x) if and only if e = and x = Xj for some j = 1, . . . , L. 
The Xi's are called the "touching points" of the two bands. We want to compute the change of 
the Chern class of each of these bands while e varies from —1 to +1. Such a change has been 
observed in numerical computations [§, || [HJ related to physical models. 

To give the result in a precise form, we need some notation. First of all, we will replace the 
N x N Hamiltonian H(e, x) by an effective 2x2 one describing the two bands under study 
in a neighbourhood O of e = 0, x = Xj. This is always possible. Using the Pauli matrices 
a = (<Ti, o~2, 03), this effective Hamiltonian can be written as 

H eff (e,x) = ^(E + (e,x) + E4e,x)) + af(e,x) , (1) 

where / takes on values in R 3 and |/| = (E + (e,x) — E_(e,x))/2. From our assumptions, it 
follows that / vanishes only at (0, Xj) in O. We then set uj = f/\f \ in the complement of the 
(0, XjYs in O. Given any surface £j in O homotopic to a sphere centered at (0, Xj) and not 
surrounding the other touching points, we set 

n(xj) = — J < uj\dQ A duo > . (2) 
This is an integer which does not depends on the choice of Uj, called the degree of the map 



uj. This index is actually related to Berry's phase. In his seminal paper |TTJ, M. Berry showed 
that if at e = the parameters x turns around Xj on a closed path 7 close enough to Xj, the 
eigenvector ip + corresponding to the +-band experiences a phase change given by 

A0 = \ ( < uj\dQ Adu > , (3) 
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where S is any oriented surface in O with boundary given by 7. For this reason, we propose 
to call n(xj) the Berry index at Xj. Our main result is the following: 

Theorem 1 The difference ACh between the Chern numbers of the subband corresponding to 
E + while e varies from —1 to +1 is given by the sum of the Berry indices at all touching points. 

In Section |3] we will compute the Berry index in some cases. The generic situation corresponds 
to / having a non vanishing Jacobian at the touching point and gives the following result: 

Theorem 2 If f has a non vanishing Jacobian Jj{xj) at the touching point (0, Xj), the Berry 
index is given by n(xj) = J^Xj)/\JAxj)\, namely it is the sign of the Jacobian at the touching 
point. 



2 Proof of the main result 

In what follows, M is an oriented smooth surface. Up to an homeomorphism, it is given by 
the Poincare construction from a polygon in the Poincare disk (see Fig. [I]). Let us assume for 
simplicity that for e 7^ 0, the spectrum of the Hamiltonian H(e, x) is simple. We then denote by 
Pi(e, x) the eigenprojection associated to the eigenvalue Ei(e, x). It is given either by \ipi >< ipi\ 
where ipi is a corresponding normalized eigenvector, or by Cauchy's formula as 
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Figure 1: The surface M 



H^) = ^-4 — % — w ( 4 ) 

lin Jr z — H (e, x) 

where T is a small circle centered at Ei(e,x) and not surrounding the rest of the spectrum. 
Since the spectrum is simple for e 7^ 0, Pi is smooth with respect to (e, x) in the complement 
of the touching points in [—1, +1] x M. For a given e its second Chern class is given by the 
closed 2-form Q 2 011 M where 

^2 = ^Tr(PidPi A dPt) . (5) 
The Chern number Ch(e) is obtained by integrating Q 2 over the surface M. That it is an 



integer is a standard result p2| , |13| , [|. Moreover, this number does not change with e as 
long as the Z-th band remains isolated from the rest of the spectrum. Finally if Pi is two- 
dimensional, namely of the form (l + au)/2 for some unit vector u in R 3 , the two-form becomes 
fi 2 = (V 2 ) < 3\dfi A dCu >. 

Now, we assume that, at e = 0, the band E + = E\ touch E- = Ei-\ at one point xq e M 
only. The general case can be obtained in a similar way. Let D be a small open neighbourhood 
of xq in M diffeomorphic to the unit disk in R 2 . Let then 7 be the oriented boundary of D 
and Mj) be the complement of D in M endowed with the same orientation as M. Therefore, 
as an oriented path 8Md = —7. For < rj < 1, we construct a new surface Md (77) by gluing 
together x (-Mo), [—77, +17] x 7 and {—77} x (Mo), where —M' denotes the manifold M' 

with opposite orientation (see Fig. ^). 

We give Md (77) the natural orientation to make it closed. As 77 — > this surface is homotopic 
to itself and becomes two copies of M73 with opposite orientations glued along 7. Moreover, for 
every 77, this surface contains no point on which the spectrum of H is degenerate. Therefore Pi 
defines a line bundle on it. Its Chern class is given by integrating Q 2 on Md(t]). Owing to the 
homotopy invariance of this Chern number, it is enough to compute it for 77 = and gives the 
obvious result 
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Figure 2: The surface Md(t]) 

L n 2 = o . 

Jm d ( v ) 

On the other hand, Mjj(t]) can be decomposed into three closed surfaces (see Fig. U), namely it 
is homologous to the sum of the following three cycles M + = {+r/} x (— M), M_ = {— rj} x (M) 
and C which is obtained by gluing together {+77} x (£)), {—rj} x (-D) and [—77, +77] x 7 with 
the corresponding orientation. Integrating Q 2 over M + gives — Ch(+1), namely the opposite of 
the Chern number of the Z-th band for e > 0. Integrating it over M_ gives Ch(— 1), namely the 
Chern number of the Z-th band for e < 0. Finally the contribution of C is precisely the Berry 
index defined by eq. (||) in the introduction, because C is homotopic to a sphere surrounding 
the touching point once. □ 



3 Computation of the Berry index 

Generic touchings - In the generic case, the function /, in eq. ([!]), expressed in a suitable 
local chart, admits a Taylor expansion about the touching point xq given by 

f(e, x + = ef + a/i + 6/2 + 0(e 2 + £ + Q , 

where £ = (£1,^2) is a small vector in R 2 and the three vectors fj are linearly independent. 
We do not change the homology class of the closed two-form Q 2 m e Q- © by modifying / 
homotopically. We can therefore ignore the remainder and deform continuously the basis fj 
to make it coincide with the canonical one up to orientation. The orientation is given by 
the sign e of < fo\}\ x f 2 >, namely the sign of the Jacobian of / at the touching point 
e = 0, x = xq. Therefore, up to homotopy, we can assume f(e,x + £) = e(e, £). Using polar 
coordinates (6 1 , 0) in R 3 , this gives Q = f/\f\ = e(cos(9), sin(6 l ) cos(0), sin(9) sin(0)), and thus 
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^2 — Esin(d)d6 A d(j). Since C is homotopic to a sphere, it follows that the Berry index in this 
case is 



n(x )=e = s g n det|^y = 

/ e=0,x=XQ ■> 



□ 



Parabolic touchings - Parabolic touching may occur non generically. Such examples have 



been encountered in the literature |9|, |10|, |bj| in connection with either quantum chaos or Bloch 
electrons in magnetic fields. By redefining the parameter e if necessary, we may assume that / 
vanishes linearly in e at x = x$. We will say that the touching is "parabolic non degenerate" 
whenever 

/(e, x + = efo + eju + £2/22 + + 0(e 2 , |£| 3 ) , 

where the three vectors {/n, ^22, /12} are linearly independent. Let e = ±1 be the orientation 
of this basis. Again by homotopy, we can neglect the remainder and replace these three vectors 
by the canonical basis of R 3 , up to the sign e. Thus / = efo + £<f(£) where 

£(0 = (&&2£i6) • 

Whenever £ varies in R 2 , g(£) describes the cone fC = {(x, y, z) £ R 3 ; x > 0, y > 0, (x + y) 2 = 
(x — y) 2 + z 2 }. This is a cone with vertex at the origin, with axis parallel to (1, 1, 0) and with 
basis given by the ellipse perpendicular to the axis. If moreover, £ describes a small circle 
7 in R 2 , <f(£) describes an ellipse similar to the basis of the cone and approaching the cone 
vertex as the radius of 7 decreases to zero. In addition, this ellipse is described twice as 7 is 
described once. Therefore we find two situations compatible with our hypothesis. Either ±/ 
belongs to the interior of this cone, and then the surface f(C) surrounds the origin twice, with 
the orientation given by e, or it does not and f(C) does not surround the origin at all. In the 
former case the Berry index is n(xo) = 2e in the latter n(xo) = 0. 

If the touching is "parabolic degenerate", namely if the three vectors {fn, fn-, /12} are 
linearly dependent, the same argument shows that the cone /C is flat so that the Berry index 



vanishes. This is the case if for instance f 12 = [14 . 
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